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Abstract. We show that if the evolution of cosmic ray sources follows that of active galactic nuclei, 
the upper limit on the diffuse flux of tau neutrinos from the Pierre Auger Observatory marginally 
constrains the proton fraction at the end of the energy spectrum. We also discuss prospects to 
uncover new physics leptophobic interactions using future Auger data. 
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The Pierre Auger Observatory is a major international effort aiming at high-statistics 
study of highest energy cosmic rays to uncover their origins and nature yj]. These cos- 
mic rays are too rare to be directly detected, but the direction, energy, and to some extent 
the nature of the primary particle can be inferred from the cascade of secondary parti- 
cles induced when the primary impinges on the upper atmosphere. These cascades, or 
air showers, have been studied in the past by measuring the nitrogen fluorescence they 
produce in the atmosphere, or by directly sampling shower particles at ground level. 
Auger is a "hybrid" detector, exploiting both of these methods by employing an array of 
water Cerenkov detectors overlooked by fluorescence telescopes; on clear, dark nights 
air showers are simultaneously observed by both types of detectors, facilitating powerful 
reconstruction methods and control of the systematic errors which have plagued cosmic 
ray experiments to date. Auger also provides a promising way of detecting ultra-high 
energy neutrinos by looking for deeply-developing, large zenith angle (> 60°) or hor- 
izontal air showers. At these large angles, hadronic showers have traversed the equiv- 
alent of several times the depth of the vertical atmosphere and their electromagnetic 
component has extinguished far away from the detector. Only very high energy core- 
produced muons survive past 2 equivalent vertical atmospheres. Therefore, the shape of 
a hadronic (background in this case) shower front is very flat and very prompt in time. 
In contrast, a neutrino shower appears pretty much as a "normal" shower. It is there- 
fore straightforward to distinguish neutrino induced events from background hadronic 
showers. Moreover, because of full flavor mixing, tau neutrinos could be as abundant 
as other species and so very low V T fluxes could be detected very efficiently by Auger 
detectors by looking at the interaction in the Earth crust of quasi horizontal v T inducing 
a horizontal cascade at the detector. 

The Pierre Auger Observatory has confirmed a strong steepening of the cosmic ray 
flux above 10 106 GeV |0, 0]. If such a steepening is due to cosmic ray interactions 
with the cosmic micorwave background, the cosmic ray arrival directions at the ener- 
gies above the suppression are expected to correlate with the nearby matter distribu- 
tion, which is quite inhomogeneous. Observation of these correlations would provide a 



doorway towards cosmic ray astronomy (provided that particle deflections in the inter- 
vening galactic and extragalactic magnetic fields are relatively small). Anisotropy clues 
have been already reported in the data: cosmic ray events above ~ 10 10 8 GeV seem 
to correlate over angular scales of less than 6° with the directions of nearby (distance 
< 100 Mpc) active galactic nuclei 0], though the evidence for such anisotropy has not 
strengthened in a more recent analysis with larger statistics In addition, upper limits 
on the diffuse neutrino flux and the ultra-high energy cosmic ray photon fraction have 
been established [5]. All these observations strongly support the hypothesis that cosmic 
ray acceleration takes place in astrophysical objects. The measurements of the variation 
of the depth of shower maximum with energy, interpreted with current hadronic inter- 
action models, favour a mixed cosmic ray (protons + nuclei) composition at energies 
above 10 8 6 GeV fl6|]. It is interesting to note, however, that uncertainties in the proton- 
air cross-section at the high energies characteristic of air showers may be large enough to 
alter this conclusion [LVD - In what follows we show that limits on neutrino fluxes provide 
complementary information to ascertain the cosmic ray nature. 

It is helpful to envision the cosmic ray engines as machines where protons are accel- 
erated and (possibly) permanently confined by the magnetic fields of the acceleration 
region. The production of neutrons and pions and subsequent decay produces neutrinos, 
gamma rays, and cosmic rays. If the neutrino-emitting source also produces high en- 
ergy cosmic rays, then pion production must be the principal agent for the high energy 
cutoff on the proton spectrum. Conversely, since the protons must undergo sufficient 
acceleration, inelastic pion production needs to be small below the cutoff energy; con- 
sequently, the plasma must be optically thin. Since the interaction time for protons is 
greatly increased over that of neutrons because of magnetic confinement, the neutrons 
escape before interacting, and on decay give rise to the observed cosmic ray flux. The 
foregoing can be summarized as four conditions on the characteristic nucleon interac- 
tion time scale T} nt ; the neutron (pion, muon) decay lifetime T n i n ^y, the characteristic 
cycle time of confinement T C ycie; the cooling time scale of charged particles T coo i; and 
the total proton confinement time T conf : (i) T int > T cyde ; (ii) x n > T cycle ; (Hi) T int < T conf ; 
(iv) tiz.pL < T CO ol- The first condition ensures that the protons attain sufficient energy. 
Conditions (i) and (ii) allow the neutrons to escape the source before decaying. Condi- 
tion (Hi) permits sufficient interaction to produce neutrons and neutrinos. Condition (iv) 
ensures that cooling processes of secondary pions, and muons have negligible influence 
on the energy distribution of the emitted neutrons and neutrinos. 

Within this working hypothesis the relative number and energy of the neutrinos and 
neutrons depend only on kinematics, which implies approximate equipartition of the 
decaying pion's energy between the neutrinos and the electron, and the relative radiation 
density in the source. On average, each interaction will produce 77 neutrinos per neutron 
with relative energy e per neutrino, i.e. 

n- {Nv) and e - {Ev) (1) 
1} -(N n ) ^ e ~W)' () 

The neutrino emissivity of flavor i is then given by: 

^J? Vi (z,E Vi ) = ^J? n (z,E n ). (2) 



Assuming flavor universality as well as £ ~ E v /E n ~ AE v /AE n and r\ ~ N a \\ v /N n we 
arrive at the neutrino source luminosity (per co-moving volume): 

SZnvizM^^Jfn&Ev/e). (3) 

In the following, we will consider a (hypothetical) source, where pion production pro- 
ceeds exclusively via resonant py^ A + scattering with fixed values rj = 3 and e = 0.07. 
Note, that the relation © derived for optically thin sources can be regarded as a lower 
limit on the neutrino luminosity as long as energy-loss processes in the source are neg- 
ligible. Next, we translate this conservative expectation into an upper limit on the extra- 
galactic proton fraction in ultra-high energy cosmic rays, exploiting the experimental 
upper limit on the diffuse flux of tau neutrinos 

For this procedure we introduce test functions of the neutron source luminosity of 
the form _S?„ test (0,£) = Sfa(E/E max) ex p(~ E/E msix ) , with an exponential energy 
cut-off E meLX that we vary between 10 8 GeV and 10 12 GeV with a logarithmic step- 
size of \og lQ E = 0.25. We adopt the usual concordance cosmology of a flat universe 
dominated by a cosmological constant with O.^ ~ 0.7, the rest being cold dark matter 
with £l m ~ 0.3. The Hubble parameter is given by H 2 (z) = Hq (Q m (l + z) 3 + £Ia), 
normalized to its value today of 70 kms -1 Mpc -1 . The time-dependence of the red- 
shift can be expressed via dz = —dt(l +z)H. The cosmological evolution of the source 
density per co-moving volume is parameterized as Sfi(z,E) = J4? (z)^Ci(0,E) , where 
the source luminosity per co-moving volume is assumed to follow that of gamma 
ray bursts [J^grb(z) = (1 +z) 4 " 8 , for z < 1; ^grb(z) = M(l + z) 3A , for 1< z < 4; 
JT grb (z) = N-iN 4 (l +z)-° A , for z > 4; with N { = 2 3J and N 4 = 5 3 2 ] and that of active 
galactic nuclei [J^agn(z) = (1 +z) 5 - , for z < 1.7; ^agn(z) = Nu, for 1.7 < z < 2.7; 

j^ GN ( z ) = N in Nfj~ z \ for z > 2.7; with N\j = 2.7 5 and N 2 j = 10 043 ] . Each neutron 
test luminosity is related to a neutrino luminosity by the ratio ©. After propagation we 
normalize the accumulated contribution of extra-galactic and cosmogenic neutrinos to 
the limit on ultra-high energy tau neutrinos from Auger. For details on the calculation 
see lUl]. The results are shown in Fig. [TJ For J^agn(z), the upper limit on the diffuse flux 
of tau neutrinos marginally constrains the proton fraction somewhat around 10 10 GeV. 

Cosmic neutrinos are also unique probes of new physics as their interactions are 
uncluttered by the strong and electromagnetic forces and can reach center-of-mass 
energies y/s ~ 245 TeV. However, rates for new physics processes are difficult to test 
since the flux of cosmic neutrinos is virtually unknown. It is possible to mitigate this 
by using multiple observables which allow one to decouple effects of the flux and cross 
section. In particular, possible deviations of the neutrino nucleon cross section due to 
new physics leptophobic interactions (e.g., production and decay of TeV-scale black 
holes) can easily be uncovered at Auger by combining information of the two different 
neutrino channels [9]. Namely, if an anomalously large quasi-horizontal deep shower 
rate is found, it may be ascribed to either an enhancement of the incoming neutrino 
flux, or an enhancement in the neutrino-nucleon cross section. However, these two 
possibilities may be distinguished by separately binning events which arrive at very 
small angles to the horizontal, the so-called tau Earth- skimming events. An enhanced 
flux will increase both quasi-horizontal and Earth- skimming event rates, whereas a large 
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FIGURE 1. Upper limits on the proton contribution in ultra-high energy cosmic rays derived from 
Auger limit on the diffuse v T -fiux in models where the cosmic ray sources evolve strongly with redshift; 
we show the case of strong cosmological evolution of the proton sources according to J%grb an d ^agn- 



black hole cross section suppresses the latter, because the hadronic decay products 
cannot escape the Earth crust. This argument requires only counting experiments and 
does not rely on measurements of shower properties. 

In summary, statistics and better experimental handles should enable us to reconstruct 
the high end of the energy spectrum, to locate the cosmic ray sources in the sky, and to 
discern the primary mass composition. Future Auger data will not only provide clues to 
the cosmic ray origin, but could enhance our understanding of fundamental physics at 
sub-fermi distances. An optimist might even imagine the discovery of microscopic black 
holes, the telltale signature of the universe's unseen dimensions. 



ACKNOWLEDGMENTS 

I would like to thank Markus Ahlers, Jonathan Feng, Haim Goldberg, Subir Sarkar, and 
Al Shapere for the most enjoyable collaborations. This work has been partially supported 
by NSF Grant No PHY-0757598 and the UWM Research Growth Initiative. 



REFERENCES 

1. J. Abraham et al. [Pierre Auger Collaboration], Nucl. lustrum. Meth. A 523, 50 (2004). 

2. J. Abraham et al. [Pierre Auger Collaboration], Phys. R ev. Lett. 1 01, 061 101 (2008). 

3. J. Abraham et al. [Pierre Auger Collaboration] , |arXivi0 906 .2189 

4. J. Abraham et al. [Pierre Auger Collaboration], Science 318, 938 (2007). 

5. J. Abraham et al. [Pierre Auger Collaboration], larXiv:0906.2347l 

6. J. Abraham et al. [Pierre Auger Collaboration], arXiv:0906.2319 

7. R. Ulrich, R. Engel, S. Muller, F. Schussler and M. Unger. larXiv:0906.30751 

8. M. Ahlers, L. A. Anchordoqui and S. Sarkar, Phys. Rev. D 79, 083009 (2009). 

9. L. A. Anchordoqui, J. L. Feng, H. Goldberg and A. D. Shapere, Phys. Rev. D 65, 124027 (2002). 



